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The construction of differential cross sections as a function of excitation energy for systems with 
a collection of low- and high-lying intrinsic vibrational modes has been attempted in the past. 
A prescription is proposed that simplifies the implementation of such calculation schemes with a 
remarkable reduction in computational time. 
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I. INTRODUCTION 

The possibility of exciting collective vibrational modes 
in heavy ion reactions has captivated the attention of re- 
searchers in the field for several decades [IJ. Of especial 
interest has been the challenge of extending the set of 
participant states (i.e. the open channels) well beyond 
the familiar low-lying surface modes and into the range 
of the nuclear giant resonances of the lowest multipolari- 
ties. Because of adiabatic considerations these processes 
required, at the bombarding energies available during the 
1970's and 80's, very short effective collision times. Only 
the sharp exponential drop of the nuclear radial formfac- 
tors was able to recreate these conditions at that time, a 
fact that was noted (and exploited) in the early studies 
of Deep Inelastic Colhsions [3, Q • 

Nowadays the accelerators provide much higher beam 
energies and thus the population of giant resonances 
can also be mediated by the Coulomb excitation mech- 
anism. Under these circumstances the long range of the 
Coulomb-coupling matrix elements forces us to incorpo- 
rate ~ in the theoretical analyses of the process - a con- 
siderably larger number of impact parameters or partial 
waves. This, naturally, increases the chances for exciting 
the high-energy part of the nuclear response. 

Historically, the multiple-phonon excitation of high- 
lying modes was actively promoted for being the source of 
characteristic structures in the experimental distribution 
of cross sections as a function of excitation energy [2]. 
The study of these patterns became, in turn, a quite 
convenient source of information for learning about the 
actual features of the giant resonances (energies of the 
modes, widths, strengths, anharmonicity, etc). 

This link motivated an important body of recent theo- 
retical work 0] where microscopic calculations for the 
structural aspects of the giant modes have been combined 



with standard reaction formalisms to yield concrete pre- 
dictions for the shape of the distributions da/dE. Unfor- 
tunately, as it was stated earlier, a much larger number of 
impact parameters are now needed to compute accurately 
the nuclear and Coulomb components of the excitation 
processes. This, together with the fact that there exists 
a distinct possibility of exciting a multiple number of col- 
lective phonons, has resulted into rather complicated and 
time-consuming coupled-channel schemes. 

Upon close inspection of the results, however, one real- 
izes that for practically the entire range of relevant par- 
tial waves the excitation probabilities are very small. It 
is also possible to conclude that, in leading order, the 
different modes can be considered as being independent 
from each other. Or, in other words, that one can ig- 
nore terms in the hamiltonian that involve simultane- 
ously the coordinates of two or more collective variables 
a\^. Notice that this is not the same as claiming that 
one works within the perturbation limit; we have already 
stressed the relevance of multi-step events. Multiphonon 
processes may occur as the excitation of the same mode 
(two, three phonons) or as the simultaneous excitation of 
two or more different modes [Ij. 

Considering the harmonic modes as being essentially 
uncoupled to each other has a significant practical ad- 
vantage; it can be exploited to design an event generator 
that allows for a much simpler, and yet accurate, method 
for constructing the differential cross sections da/dE. We 
shall describe this idea in detail in the following sections. 

Such an approximation scheme is of course bound to 
fail for the very central impact parameters. However, 
these partial waves are unlikely to participate directly in 
the population of the inelastic channels explicitly taken 
into account. Within the semiclassical formalism their 
contribution is, in fact, strongly suppressed by the global 
absorption associated with the imaginary part of the op- 
tical potential. To take into account properly the grazing 
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impact parameters may perhaps require a cornplete pro- 
cedure as the one exploited by the authors of @. Notice 
that this would be necessary, at worst, for only a narrow 
window of impact parameters; for most of the relevant 
range (extending up to hundreds of fermis) the method 
we are proposing in this paper is appropriate and it is 
precisely here that the major savings of computational 
time can be achieved. 

We elaborate further on these considerations in 
Sect. II, where we also recall a previous work that proves 
quite essential to the development of our proposal. In 
Sect. Ill we explain what is truly specific about the tech- 
nique for generating the folding of probabilities that we 
need. In Sect. IV we take, as an example, the reaction 
""^Ar -I- 208p|^ per nucleon and include two 

low-lying and three high- lying resonances. The purpose 
is here mainly to show that the function da / dE{E) in- 
deed reflects correctly the assumed input to the problem. 
We reserve Sect. V for a brief summary of the contribu- 
tion and some closing remarks. 

II. BRIEF BACKGROUND 

A few years back we developed a general formalism for 
the excitation of a single collective vibrational mode 0] . 
We begin the presentation mentioning this reference be- 
cause it proved quite essential to motivate the present 
contribution. Exploiting well-tested approximations we 
were able in [Ti] to propose a semiclassical prescription to 
estimate the dynamical effects associated with a spread- 
ing width r of the mode. Modifications in the predicted 
cross sections due to the presence of an eventual anhar- 
monicity, manifested by an apparent ratio of state ener- 
gies e 

e(2 phonon state) ^ 
£(1 phonon state) ' 

were also investigated. We should mention here that the 
development of this program - aiming mostly to the de- 
scription of single- and double-phonon giant resonances - 
was done within the framework of perturbation theory. 

Without entering into details in the implementation of 
ref. let us briefly recall what was the input of the cal- 
culation scheme and what could be obtained as a result. 
Given a single collective vibrational mode of multipolar- 
ity A, energy hwx, width T\ and anharmonicity v\ a run 
of the program densely sampling impact parameters p 
belonging to the interval [pmim Pmax] generated the to- 
tal differential cross section da / dE\{E). The energies of 
the transitions — > 1 phonon I — > 2 phonon, affected as 
they are by the anharmonicity factors and the spreading 
widths, introduced interesting dynamical consequences 
which were the main object of investigation in Q. 

To construct the differential cross sections da/dE (E) 
it was necessary to define a procedure to distribute the 
probabilities of inelastic transitions for the one- and two- 
phonon vibrational states over the relevant excitation en- 



Ar + Pb @ 40 MeV/nucleon 
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FIG. 1: (Color online) Treatment of the absorption. The fig- 
ure on the top shows the behaviour of T(b) as function of 
the impact parameter b. In the middle frame, the excitation 
probability of a single phonon is plotted as function of b for 
the three states reported in the legend. The functions ob- 
tained as a product of the two previous quantities is plotted 
in the figure at the bottom. 

ergy ranges. This prescription is quite analogous (except 
for its generalization from one mode to several modes) 
to the one we later use for our event generator. We can 
thus defer its presentation to the next Section. 

The effect of the absorption was taken into account 
in [rj by means of a multiplicative depletion factor that 
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rapidly falls to zero as the overlap between the reacting 
nuclei increases for the lower impact parameters. This 
is - following standard practice - constructed from an 
integral along the trajectory r{t), 

T{b) = exp{--J W{r{t')dt'} (2) 

of the imaginary part of the optical potential, W (see for 
instance, j8| and refs. quoted therein). For the current 
application the absorptive component was chosen follow- 
ing the prescription of refs. 0jQ. 

The projection into the subspace of explicitly consid- 
ered channels is now able to reverse the tendency of the 
second-order amplitudes to yield too large probabilities 
for impact parameters at or inside the grazing distance. 
Notice that this uncomfortable situation would not be- 
come obvious when solving for coupled-channel ampli- 
tudes (even if the numerics may be equally absurd) be- 
cause of a prescribed conservation of the norm by the 
integration algorithm. 

The interplay between these contrasting effects can be 
appreciated in Fig. 1 for the reaction '^^Ar -I- ^°*Pb at 
40 MeV per nucleon. The figure shows the attenuation 
factor that defines what fraction of the contribution of 
a given impact parameter is actually retained in the in- 
elastic channels. The top frame shows the transmission 
coefficient T(b) as a function of impact parameter and 
puts in evidence that the transition from T{b) « 1 to 
T{b) w indeed occurs over an interval Ah that spans 
only a couple of fermis. The middle frame displays the ex- 
citation probabilities as a function of the impact param- 
eter in three situations. These are, namely, the low-lying 
quadrupole mode (full curve), the giant quadrupole res- 
onance (GQR, dashed curve) and a high energy octupole 
resonance (HEOR, dotted curve). The lower frame cor- 
responds to the actual distribution da/db, constructed 
from the information displayed in the two plots above. 

III. FORMALISM 

We consider a pair of reacting heavy ions that accumu- 
late a number N of intrinsic surface vibrational modes. 
Each one is characterized by their excitation energy fkUi , 
multipolarity A.;, strength Pi and a width F^. That is, 

hu,, A,, A, ,r, with l<i<N. (3) 

We proceed immediately to define two separate groups 
of these modes: high- and low-lying modes. The mean 
reason for establishing the subdivision has to do with the 
way people in the field have traditionally dealt with their 
widths. The differentiation is actually an old story that 
dates back to the 70 's and the use of surface vibrational 
models to describe specific dynamical features of Deep 
Inelastic Reactions [1, It aims to refiect two clear 
experimental facts: 



• Low-lying modes (excitation energies Hux < 5- 
6 MeV) are such that at the zero- and one-phonon 
level are sharp and display no width. Clearly the 
energy range quoted above is only qualitative. The 
modes we have in mind are those, for example, 
known - in the harmonic oscillator terminology for 
even multipolarities - as AiV=0 excitations. At 
the two-phonon level they show a spread which is 
mainly associated with the anharmonicity of the 
mode. Let us be a bit more explicit. Suppose we 
have a low-lying quadrupole mode with hjj\^2 ~ 4 
MeV. At about double that excitation energy it is 
found a multiplet of states 0+, 2+, 4+ spanning an 
interval of energy that we shall call AEx=2 ■ Typ- 
ically this quantity has an order of magnitude of 
about 0.5 MeV. 

Finally, the last piece of experimental evidence that 
one wishes to incorporate in the formalism is that 
there are practically no known three-phonon states 
associated with low-lying modes. This can be for- 
mally done by assuming that, at the three-phonon 
level, the mode assumes a width that equals the 
separation energy hiox. 

All of these features are best implemented by as- 
cribing to the mode an energy-dependent width 
(see below). 

• High- lying modes (i.e. giant resonances). In this 
situation the zero-phonon state is taken to be 
sharp, while at the one-phonon level the state dis- 
plays the well-known spreading width F^. The dis- 
tribution of the excitation amplitudes to higher lev- 
els is achieved by a straightforward folding (see be- 
low). 

The prescription that emerges from the two items 
listed above may appear, at first, difficult to grasp. How- 
ever it has led to practical conclusions in the treatment of 
Deep Inelastic Collisions that are in very good agreement 
with the experimental evidence. Obviously one could 
come up with different but somewhat equivalent oper- 
ating procedures. It is easier and reasonable, however, to 
adhere to this established practice since the details of its 
implementation have already been described and tested 
in the literature (see for example ref.jl]). A practical re- 
minder of the spreading prescription as it is applied for 
low- and high- lying modes is summarized in Table 1. 

The event-generator works in practice just like any 
other similar device in a wide variety of physics sub- 
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TABLE I: Width prescription as it is used in the construction 
of the cross section shown in figures 2 and 3. 
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FIG. 2: (Color online) Width prescription as applied to low- 
and high-lying modes. The distributions of the total differ- 
ential cross section da / dE{E) is constructed as an example 
for two separate quadrupole modes with fi£j=3 and 12 MeV. 
In this figure we can appreciate the different prescriptions for 
low- and a high-lying modes summarized in Table I. The giant 
quadrupole state (r(n = 1)= 5 MeV) clearly shows the sec- 
ond phonon with a wide and increasing width. The low-lying 
mode (Ai?= 1 MeV) displays, on the other hand, a much 
narrower structure at the two-phonon level and has the three 
phonon structure practically washed out (not visible in the 
figure). A large number of events and an impact parameter 
range of 100 fm was used for this illustration. 



jects. We take a given impact parameter and consider 
a very large number of possible "events", N^v In every 
instance one generates the probability of occurrence of 
each independent mode by "throwing a Poisson dice" N 
times, in agreement with an average number of phonons 
{rii). This means that the random number generator is 
designed to return a number of phonons for each of 
the independent modes consistent with the law 



P{.ni) 



exp{-(ni)} . 



(4) 



The differential cross section we search for is to be con- 
structed for a large range of relevant impact parameters 
[pii P2\i that are sampled with a uniform interval Ap. 
The average number of phonons for the i^^ independent 
collective surface mode is previously calculated, for each 
impact parameter p^, using the formalism developed in 
ref. These figures are collected in an auxiliary data- 
file that has the structure 



Pi, ),..., (rii^jv) 



Pk, (ni=i ),.-., (?^i=jv) 



P2, (nj=i)...., (ni=iv) • 



(5) 



The smaller impact parameter, pi should be such that 
any possible contribution from it is definitely eliminated 



by the absorption. Similarly, one should verify that con- 
tributions from p > P2 can also be neglected. 

Adding to this prepared data-set the characteristic in- 
formation that specifies the different vibrational modes 
i one is ready to run the event generator and construct 
da{E)/dE. 

For each impact parameter in the file ([5]) the probabil- 
ity assigned to the current event is, naturally. 



N 



(6) 



It is in the energy scale that we have to be careful with 
the character of low-lying vibrational state or giant res- 
onance of the particular mode i. Following the prescrip- 
tion summarized in Table 1 we assign a spread F^ to the 
contribution to the excitation energy, e^, of this mode. 
This quantity is thus defined as 



Hi huji + Q{0, Si) , 



(7) 



where 0(0, T,i) is a random number obtained from a nor- 
mal distribution with zero centroid and standard devia- 
tion Ei w Fi/2.3. Slightly different prescriptions could 
once more be obtained by replacing Q hy a similar type 
of random number generator, but these choices are not of 
much consequence at the level of approximation we have 
chosen to maintain. 

The total excitation energy for the collection of N in- 
dependent modes is then simply given by 



E 



N 

E 

1=1 



(8) 



absissa to which - in a properly designed histogram - we 
accumulate a weighted version of the probability previ- 
ously given in eq. ^ . The proper units for the differen- 
tial cross section as a function of energy are obtained 
by multiplying that number by {20TTpkAp)/{AENev), 
where pk is the impact parameter sampled at this mo- 
ment, and Ap, AE are the extents of the impact param- 
eter mesh and the energy mesh respectively. The final 
result is then given in mb/MeV and a drawing of this 
histogram represents the predicted distribution of cross 
section da{E)/dE. 

To test the validity of our method, a comparison with 
a more sophisticated method like the one of ref. Q is in 
order. In this approach one starts with a Hartree-Fock 
plus Random Phase Approximation calculation in order 
to identify the most collective one-phonon states. For 
each of these chosen states one calculates their transition 
densities and the corresponding form factors. These are 
used in a semiclassical coupled channel scheme to deter- 
mine the excitation probabilities for all the possible one-, 
two- and three-phonon states that one can construct. 

In Table II we compare some results for the two ap- 
proaches. The calculations are performed for the system 
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^°Ga. + 208pb at 50 MeV per nucleon. We take a simple 
example where only two one-phonon states are consid- 
ered as input: the low-lying 3~ state (E=4.9 MeV) and 
the GQR (E=16.9 MeV). The range of impact parameter 
used in the calculations (15-100 fm) corresponds to the 
peripheral region where the Coulomb interaction yields 
the most important contribution. This is also the re- 
gion where the excitation probability distribution are of 
a Poisson type. The results of the two methods are very 
close. 

So, one can envisage a calculation performed in two 
steps: Make use of the method of ref. Q in the inner 
region where the nuclear interaction plays an important 
role and then use our novel approach in the peripheral 
region for large impact parameters which is the most time 
consuming part. 

Finally, we would like to stress that in the case one 
wants to take into account the contribution of anhar- 
monicities, our method does not apply. 

IV. APPLICATION 

We now proceed to illustrate the possibilities of the event 
generator with an application to the specific reaction 
38^j. _|_ 208p|3, j^i- 4Q ]y[gY pgj. nucleon. We take two low- 
lying modes (one quadrupole, one octupole) with ener- 
gies huj=3 and 5 MeV respectively, and a common value 
AE—1 MeV. The deformation parameter assumes for the 
low-lying modes a value (3=0.1 . Three giant resonances 
are then added; a dipole, a quadrupole and an octupole 
mode with energies hu!=18, 17, and 31 MeV and widths 
r=6, 6, and 8 MeV respectively. The effective spread of 
all these modes is found in all cases following the prescrip- 
tions described in Table 1. Accumulating the cross sec- 
tion for a range of impact parameters from pmi„=12 fm 
to Pmax=100 fm in steps of 0.5 fm we obtain the distri- 
bution shown in Fig. 3. 

The characteristics of the curve refiect the assumptions 
made and, in addition, are reminiscent of what is ob- 
tained by the time-demanding method of ref. [5|] in sim- 
ilar circumstances. Just to stress the advantage of the 
proposed method let us note that one can save about 
two orders of magnitude in computing time by using the 
event generator. 
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0.14 xlO"^ 
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TABLE II: Total excitation cross section (in mb) calculated 
with the events generator method (second column) and with 
the method of ref.[^ (third column) for the states shown in 
the first column. 



V. CONCLUDING REMARKS 

The developments described and the results shown in this 
Report were motivated by sheer necessity. In fact, while 
the complete calculations performed by Lanza et al. are 
very important the absolute times involved in the com- 
putation of cross sections by the procedure described in 
ref. Q are quite large. This becomes even more critical 
when one takes into account that the practical imple- 
mentation of such a prescription involves a considerable 
amount of leeway that can only be sorted out by trying 
alternative calculations with equally acceptable sets of 
parameters. That appears to be the only sensible way 
to learn, by gaining experience, how the various input 
numbers do indeed affect the calculated distributions of 
cross section. Even if a final presentation is contemplated 
with the full procedure of [5] a number of previous cal- 
culations exploiting the event generator would be - no 
doubt - very convenient to prepare the ground. 

We consider that the use of an event generator like the 
one described in these pages is highly advisable for that 
class of problems that can use suggestive results to judge 
the soundness of the answers they provide. Actually this 
should be done even before embarking in sophisticated 
schemes without the proper benefit of an educated intu- 
ition. 
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FIG. 3: Constructed distributions of the total differential 
cross section da/dE{E) for the reaction ^*Ar -I- ^"^Pb at 40 
MeV per nucleon. Some 10* events per impact parameter 
have been generated and their contributions accumulated to 
construct this figure. The calculation time is however a few 
minutes CPU. The example is described in detail and the 
relevant parameters given in the text of Sect. IV. 
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